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Abstract. We present our description of neutrino induced charged current quasielastic scattering (CCQE) in nuclei at energies 
relevant for the MiniBooNE experiment. In our framework, the nucleons, with initial momentum distributions according to the 
Locai Fermi Gas model, move in a density- and momentum-dependent mean field potential. The broadening of the outgoing 
nucleons due to nucleon-nucleon interactions is taken into account by spectral functions. Long range (RPA) correlations 
renormalizing the electroweak strength in the medium are also incorporated. The background from resonance excitation events 
that do not lead to pions in the final state is also predicted by propagating the outgoing hadrons with the Giessen semiclassical 
BUU model in coupled channels (GiBUU). We achieve a good description of the shape of the CCQE distribution extracted 
from data by MiniBooNE, thanks to the inclusion of RPA correlations, but underestimate the integrated cross section when 
the standard value of = 1 GeV is used. Possible reasons for this mismatch are discussed. 
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INTRODUCTION 

By charged current quasielastic scattering (CCQE) one usually understands the reaction in which the elementary 
process 

Vi{p)+n{k)^l-{k')+p{p') (1) 

takes place inside the nucleus. CCQE, the largest reaction channel for Ey <2 GeV, is of cardinal importance for 
oscillation experiments that rely on the detection of muons (electrons) in disappearance (Vg appearance) searches. 
It is also the channel that can be more reliably used for a kinematical neutrino energy reconstruction, an indispensable 
exercise for a precise determination of oscillation parameters in long-baseline accelerator experiments. Moreover, 
quasielastic scattering is interesting by itself and has, indeed, been carefully investigated with electron beams both 
experimentally and theoretically with a large variety of models IH Q, S m- With neutrinos it is possible to study 
different properties of the nuclear response in the axial sector that are not (easily) accessible in electron scattering 
experiments. Provided that nuclear effects are under control, CCQE could be a source of information about the nucleon 
axial form factor, often parametrized as 

Fa{Q')=8a(i + ^J , (2) 

where — —{k — k')'^ and Ma is the so called axial mass. 

The definition of CCQE given above already implies the assumption that the neutrino-nucleus interaction takes 
place predominantly on a single nucleon (impulse approximation) although at small momentum transfer |q| = |k — k'| 
(in the Laboratory frame), collective effects involving several nucleons should play a role. On the other side, as the 
excitation energy ((0 = ko — k^^) increases, inelastic channels 

Viip)+Nik)-^1-Ìk')+X{p') (3) 

with X = (ph)A^, 7zN, ... start to open. As these processes are not always identified experimentally (pions can 
be absorbed and the nuclear products may not be ali detected), they cannot be separated from CCQE in a model 
independent way. 

The MiniBooNE experiment, running with (Ey) ^ 750 MeV on a CH2 target, has collected the largest sample 
available so far for low energy CCQE [5]. After subtracting the non CCQE background, mainly from A(1232) 
excitation, using the NUANCE event generator 1^, the CCQE data set was analyzed with the relativistic Global Fermi 



Gas model of Smith and Moniz (SM) fì\\. The shape of the muon angular and energy distributions averaged over the 
Vfi flux {d(j/dcos9fidEn) could be described with rather standard values of the Fermi momentum pF = 220 MeV 
and binding energy Eb = 34 MeV, but restricting the available phase space for the final proton by means of an ad hoc 

parameter fc = 1.019 ±0.011 such that p'^-^ = k{^M^ + pj- ~ (0 + Eb), and taking Ma = 1.23 ±0.20 GeV [5]; this 

value ofMA is in agreement with the K2K result = 1.2±0.12GeV iSj] on '^O, but considerably higher than the one 
obtained from v^-deuterium data = 1 .01 37 ± 0.0264 GeV lH or by NOMAD at high energies (3-100 GeV) also on 
'^C [Ma — 1 .05 ± 0.02{stat) ± 0.06{syst) GeV] fl^. A recent reanalysis using charged current single pion production 
(CCl;r) data to adjust the Monte Carlo simulation employed to subtract the background obtains K = 1.007 ± 0.007 
and Ma — 1.35 ±0.17 GeV lITlll . While the first shape-only fit falls short compared to the measured integrated cross 
section [5], the second one underestimates it only by 10 % ifTll] . 

While such a modified SM model might be convenient to parametrize the CCQE cross section using a small number 
of parameters, it is important to understand the MiniBooNE CCQE data with more realistic nuclear models that 
implement the knowledge gathered through many years of research in electron-nucleus scattering. Besides, the fact that 
the description changes with the background subtraction procedure indicates the need of making theoretical predictions 
that could be compared to (more) inclusive and less model dependent data. 



THE MODEL 

The scattering amplitude for the elementary process [Eq. ([T]|] is proportional to the product of the leptonic and hadronic 
currents 

Jé= r- ' (4) 
V2 

While the charged-current leptonic current is given by the Standard Model, the hadronic one can be written in terms 
of form factors that contain the Information about nucleon properties 



Ja = u{p') 



u{p). (5) 



The vector form factors are obtained from electron scattering |9]; using PCAC Fp can be expressed as a function 
of Fa, given by Eq. ^ with gA~^ .267; for Ma we adopt a value of 1 GeV, consistent with the world data. 

Our description of the CCQE reaction on nuclei is based on a Locai Fermi Gas model, i.e. at each space point 
the initial-nucleon momentum distribution is given by a Fermi sphere /(r,p) = &{pF{r) — |p|) with radius pF{r) = 
[j7r^p(r)]'/^, where p(r) is the empirical nuclear density. A Pauli blocking factor for the final nucleon fi>auii = 
1 — 0(/7f (r) — IpI) also applies. Such a simple model already incorporates a space-momentum correlation which 
is absent for the Global Fermi Gas |12], and provides a framework where more elaborated many body dynamics 
can be naturally incorporated. In contrast to the Constant binding of the SM model, here ali the nucleons, initial and 
final, are embedded in a density and momentum dependent potential V(p,r) whose parameters have been fixed by 



ityand m 

proton-nucleus scattering data lll2[[l3ll . As a consequence, the nucleons acquire effective masses meff(p,r) given by 
p2 + «2, + y (p, r) = .y/p2+m2g:. 
The presence of nucleon-nucleon (NN) interactions inside nuclei implies that nucleon propagators are dressed with 
complex selfenergies E. This leads to spectral functions 

S(„)^_L ImE(/7) 

n [p^ -mi- ReZip)]'- + [lmL{p)? ' 



As most of the nucleons in the nucleus can be described as occupying single-particle states in a mean field potential 
we can neglect NN interactions for the initial nucleons (holes) and take ImE = 0. Then, Sij{p) 5{p^ — ni^^) and we 
recover the description of the initial state outlined in the previous paragraph. On the contrary, for the final nucleons 
(particle states) NN interactions should be considered. For this purpose in Eq. (|6]l we take ImE = —^/{jP^Yì^owÌPt''), 
with the collisional broadening Fcoii — POnnVysi fixed according to the parametrizations of the Giessen Boltzmann- 
Uehling-Uhlenbeck (GiBUU) framework ifTsIl . As for ReE, it is obtained from ImE with a once-subtracted dispersion 

relation demanding that at the pole p^^^^^ = \/ + '«^g. More details can be found in Ref. Il 



To complete our model we take into account that inside nuclei, the strength of the electroweak couplings may change 
from their free nucleon values due to the presence of strongly interacting nucleons iflóll . In the nuclear medium, the 
axial couplinggA is renormalized in the same way as the electric field of a dipole is screened in a dielectric medium UtIi 



gA 



(7) 



where Xd is the elementary dipole susceptibility and g' the Lorentz-Lorenz factor whose classical value is g' ~ 1/3. 
This quenching of gA in nuclear Gamow-Teller j3 decay is well established experimentally: {gA)eìi/gA ~ 0.9 lUSll 
and was first applied to CCQE scattering by Singh and Oset iflóll . Such medium polarization effects involve several 
nucleons in the nucleus and are therefore important at low |q| where the space resolution of the probe is large compared 
with the average NN separation. This corresponds to the region where MiniBooNE data exhibit a reduction with respect 
to the SM model handled in the analysis by introducing the K parameter Following Nieves et al. [|19l] . we modify the 
lepton-nucleon interaction by an infinite sum of particle-hole (ph) states (RPA), as illustrated in Fig. [T] which interact 
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FIGURE 1. Long range correlations in CCQE scattering. Solid lines pointing to the right (left) denote particle (hole) states. The 
doublé line stands for the A(1232). 

with an effective potential cast as 

V = TxT2(y\aÌ[qiqjVL{q) + [dij - qiqj)VT{q)]+ c^f X^Xi; 
Vl (Vt) contain explicit 7t (p) exchange 



(8) 



filNTl ì I ^^71 '"JT 



Ai - mi ^ 



-8 } ,Vt- { Cp I ^2 _ 2 



(9) 



and a short range part effectively included in the phenomenological Constant g' with values in the range ^' = 0.6 ± 0. 1 . 
Details about couplings and cutoff parameters A^ p can be found in Ref. Ilioll . No meson exchange is directly associated 
with the scalar term in Eq. (O, assumed to be density dependent 



P(0) 



•/(in) 



P(0) 



■l{ex) 



(10) 



where the parameters /'("') = 0.33, /'f''-'^) = 0.45 (and cq = 380 MeV fm^) are tuned to describe collective nuclear 
excitations ||20ll . The RPA sum also includes A-hole excitations as shown in Fig.[T] The ph-Ah and Ah-Ah interactions 
can be obtained by replacing C (T) with the spin (isospin) 1/2 ^ 3/2 transition operators S (T) in Eq. ^ and fyN-a 
by fàNn in Eq. Explicit expressions for the RPA corrections to the hadronic tensor are given in Appendix A of 
Ref. ìli. 

This RPA approach, built up with the single-particle states of the Locai Fermi Gas is simpler than other more 
sophisticated methods such as continuum RPA and applies only to inclusive processes, but it incorporates explicitly n, 
p exchange and A-hole states and can be naturally inserted in a unified framework to study different neutrino-induced 
reactions like inclusive quasielastic scattering, nucleon knockout and pion production |21]. Moreover, it had been 



Only the terms that contiibute to CCQE are shown. 



successfuUy applied to photo- and electro-nuclear reactions ll22l [U and allows to describe simultaneously inclusive 
muon capture on '^C and the low energy LSND CCQE measurements (with Ve and v^) [19]. 

Finally, in order to compare to model-independent data it is necessary to include the contributions from the processes 
in Eq. that look like CCQE events in the detector. The main source of such a background is pion production from 
A(1232) excitation (v^A^ M A) followed by absorption (AA' NN). Pion final state interactions in the nuclear 
medium are treated with a semiclassical BUU model in coupled channels (GiBUU, see Refs. ll23ll24ll for details). 



RESULTS 



In Fig. 12] we present the predictions of our CCQE model on '^C averaged over the MiniBooNE flux fl^ for 
four different distributions. It is important to stress that, while the upper two correspond to directly measurable 
quantities (energy and scattering angle of the outgoing muons), the lower ones can only be experimentally obtained 
by reconstruction, assuming that the target nucleon is at rest (3]. Further corrections can be made by mapping 
reconstructed to true energy with the help of a reaction model, as done by MiniBooNE with their Monte Carlo 
simulation ll26ll . The plots reveal that the nuclear many-body corrections taken into account reduce the cross section. 
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FIGURE 2. Differential cross sections for the CCQE reaction (TJ on '^C and with l = pi averaged over the MiniBooNE flux as a 
function of the cosine of the outgoing muon angle (upper left), its energy (upper right), the (true) four-momentum transfer squared 
(lower left) and the (true) neutrino energy (lower right). Dotted lines represent the Locai Fermi Gas model with Fermi motion and 
Pauli blocking. In the dash-dotted lines the nucleons are exposed to the mean field potential while the dashed ones also incorporate 
spectral functions for the outgoing nucleons. The full model with long range (RPA) correlations is denoted by solid lines. 

The effect of the spectral functions is rather small for these observables but the long range correlations cause a 
considerable reduction at forward angles (cos 9^ < 0.8) and low < 0.3 GeV^. 

As discussed above, model-independent comparisons to data must include the fake CCQE background. Our predic- 
tion for this background from A excitation and the full CCQE-like yield is shown in Fig. [3] for dcr/dQ^. In total we 



obtain a fraction of fake CCQE over the total CCQE-like events of 10 %, slightly smaUer than the prediction of the 
MiniBooNE Monte Carlo: 12 % (9.4 % CCl7r+ resonantplus 2.5 % CCItt") Ili]. 
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FIGURE 3. distribution of CCQE-like events averaged over the MiniBooNE flux (solid line). It is given by the sum of the 
pure CCQE contribution (dashed line) plus fake CCQE events where the A is excited but no pion is produced (dash-dotted line) 

CCQE-like data from MiniBooNE are not yet available. Nevertheless, it is useful to compare our pure CCQE 
predictions with the results of the MiniBooNE analysis that describe the shape of {da /dQ^) with a modified SM 
model iHdll]. Such a shape-only comparison for the distribution is presented in the left panel of Fig.|4] The two 
sets of (k-,Ma) values, thenew K= 1.007, = 1.35 GeV [denoted as (1)] and the originai 1.019,Ma 1.23 GeV 
[denoted as (2)] are shown, both actually leading to very similar shapes. The comparison reveals that the RPA 
correlations, which reduce the size of the peak with respect to the tali, bring the shape of our distribution dose to 
those extracted from data while keeping Ma = 1 GeV. On the other side, as it is clear from Fig.|2l the integrated CCQE 
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FIGURE 4. Shape of the distribution for the v^-induced CCQE reaction on '^C averaged over the MiniBooNE flux. On the 
left panel, the prediction of our model without (solid line) and with RPA correlations (dashed line) is compared to the modified SM 
model with k = 1.007, Ma = 1.35 GeV (1) and K = 1.019, Ma = 1.23 GeV (2). On the righi panel, the modified SM model (2) 
is confronted with the present model (including RPA) evaluated for different values of Ma ■ Ali curves are normalized to the same 
area. 

cross section obtained with our full model ({a) = 3.2 x 10^^^*^ cm") is smaller than the one obtained within a standard 
Fermi gas model with Ma = 1 GeV. Instead, MiniBooNE gets a considerably larger value: (a) — 5.65 x 10^^^ cm^ 
with an eiTor of 10.8 % |11]. We have explored the sensitivity of the results to some of the uncertain magnitudes in 
the model. In particular, changing g' , whose contribution to the correlations is the largest at low Q^, within acceptable 
values g' = 0.5 — 0.7 leaves the shape of the distribution practically unchanged; the impact on the integrated cross 
section is also small: {a) = 3.1 — 3.4 x lO^-'^ cm^. Increasing the value of Ma causes an increase in the total cross 
section, but at the same time, the description of the shape gets worse, as illustrated on the right panel of Fig.|4] 



CONCLUDING REMARKS 



The theoretical model for the CCQE reaction presented here incorporates important many body corrections to the basic 
Locai Fermi Gas picture like spectral functions and long range RPA correlations. The latter are important at low 
where coUective effects play a role. We find that a good agreement with the shape of the CCQE distribution 
extracted from MiniBooNE data with Ma — 1 GeV, which is favored by early neutrino data, by the analysis of 
pion electroproduction dose to threshold and by a recent neutrino experiment at high energies (NOMAD). However 
our description clearly underestimates the MiniBooNE integrated CCQE cross section. The situation is common to 
other models that take into account short-range correlations |27] or apply the phenomenological scaling function 
extracted from electron scattering [28]. Other many-body mechanisms like meson exchange currents might add some 
additional strength and should be investigated. One should also bear in mind that the MiniBooNE result is not at ali 
model independent. It relies on simulations to determine the neutrino beam and to subtract the fake contributions 
to the CCQE-like cross section. The plot of the MiniBooNE integrated CCQE cross section as function of the 
reconstructed neutrino energy together with NOMAD data (Fig. 6 of Ref. [11]) shows that, to make both experimental 
results compatible, the CCQE cross section would have to exhibit an unusual behavior, decreasing 20 — 30 % in the 
=2 — 4 GeV region to saturate afterwards. Further joint theoretical and experimental work is necessary reconcile 
available CCQE theoretical calculations and experimental values. 
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